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ABSTRACT 

We have made radio continuum, H i and X-ray observations in the direction of the 
unidentified EGRET source 3EG J1410—6147, using the Australia Telescope Gompact 
Array and the Ghandra X-ray Observatory. The observations encompass the supernova 
remnant (SNR) G312.4-0.4 and the two young pulsars PSRs J1412-6145 and J1413- 
6141. 

We derive a lower distance limit of 6 kpc to the SNR, although interpretation of 
positive velocity features in the H i spectrum may imply the SNR is more distant than 
14 kpc. PSR J1412-6145, with an age of 50 kyr, is the pulsar most likely associated with 
SNR G312.4-0.4. X-rays are not detected from either pulsar and diffuse X-ray emission 
near the bright western edge of the SNR is weak. Although there is circumstantial 
evidence that this western region is a pulsar wind nebula (PWN), the embedded 
pulsar PSR J1412-6145 is apparently not sufficiently powerful to explain the radio 
enhancement. The origin of the electron acceleration in this region and of the y-rays 
remain unidentified, unless the distance to PSR J1413-6141 is at least a factor of 3 
lower than its dispersion measure distance. 
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1 INTRODUCTION 

Since the discovery of the first unresolved y-ray source some 
30 years ago, a large number of high-energy sources re¬ 
main without obvious counterparts at other wavelengths. In 
the 3rd EGRET catalogue of y-ray sources (Hartman et al. 
1999) there are ~170 sources which are as yet unidentified. 
There have been a variety of approaches used in attempting 
to determine the nature of these unidentified sources. How¬ 
ever, due to very large positional uncertainties for EGRET 
sources (of order 1°), identifying counterparts at other wave¬ 
lengths is problematic, particularly in the Galactic Plane 
where the surface density of sources is high. The most pow¬ 
erful technique for identifying sources at other wavelengths 
is a detailed multiwavelength study of individual sources 
(e.g. Roberts et al. 2001). 

There are two y-ray sources in the southern Galac¬ 
tic Plane near longitude 312°. Observations of 2EGS 
J1418—6049 at radio and X-ray wavelengths have shown 
that it is likely associated with the young radio and X- 
ray pulsar PSR J1420-6048 (Roberts, Romani & Johnston 


2001). A hard X-ray source embedded in a radio nebula, 
possibly a pulsar wind nebula (PWN) may also contribute 
to the y-ray flux (Roberts et al. 1999). The second y-ray 
source, 3EG J1410—6147 (2EGS J1412—6211), was tenta¬ 
tively linked with the supernova remnant (SNR) G312.4—0.4 
as the most plausible radio counterpart (Hartman et al. 
1999). This association was discussed by Case & Bhat- 
tacharya (1999) but they did not reach any definitive con¬ 
clusions, merely speculating that the y-ray emission could 
be due to either emission from a PWN within the SNR or 
the acceleration of cosmic rays by the supernova shock. 

SNR G312.4—0.4 was discovered by Caswell & Barnes 
(1985) who estimated a lower distance limit of 3.8 kpc from 
H 2 CO absorption measurements. SNR G312.4—0.4 appears 
in the Molonglo Observatory Synthesis Telescope (MOST) 
SNR catalogue (Whiteoak & Green 1996) and in a de¬ 
tailed study of this region with the MOST at 843 MHz by 
Whiteoak, Gram & Large (1994). At this frequency, the SNR 
has a horseshoe-like morphology, with weaker emission in 
the south. There is a faint arc of emission about 12' to the 
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Name 

RA 

J2000 

(hms) 

Dec 

J2000 

(dms) 

P 

(s) 

P 

(xlO-i®) 

(s s-l) 

Age 

(kyr) 

DM 

(pc cm~^) 

d 

(kpc) 

E 

(xlO®"^) 
(erg s-l) 

PSR J1407-6153 

14 07 57 

-61 53 59 

0.7016 

8.85 

1250 

645 

9.7 

0.12 

PSR J1412-6145 

14 12 08 

-61 45 29 

0.3152 

98.65 

50.6 

515 

7.8 

12.4 

PSR J1413-6141 

14 13 10 

-61 41 13 

0.2856 

333.44 

13.6 

677 

10.1 

56.5 


Table 1. Properties of three radio pulsars in the error box of 3EG J1410—6147 


south of the remnant, presumably a blow-out of part of the 
shell. Observations have also been made at 4.5 and 8.55 GHz 
with the Parkes radio telescope of both total and linearly 
polarised intensity (Whiteoak 1993). The eastern boundary 
and the western enhancement of the SNR were found to be 
strongly polarised at 8.55 GHz (~10—20%), while lower lev¬ 
els of polarisation were found at 4.5 GHz (~2—5%). This is 
potentially indicative of a PWN. Figure 1 shows the MOST 
image of the SNR overlaid with the 68%, 95% and 99% con¬ 
fidence contours for the position of 3EG J1410—6147. 

There are three pulsars within the 99% confidence limits 
for the position of 3EG J1410—6147, all recently discovered 
by the Parkes Multibeam Survey (Manchester et al. 2001 
and http://www.atnf.csiro.au/research/pulsar/catalogue). 
Table 1 lists their position, period (P), period derivative 
(P), characteristic age (given by P/2P), dispersion mea¬ 
sure (DM), distance (d) determined from the DM using the 
electron density model of Gordes & Lazio (2002), and spin- 
down energy (P). The DM determined distances are known 
to be problematic in this part of the Galactic plane (John¬ 
ston et al. 2001) and are likely to be uncertain by a factor 
of 2. The two young pulsars both lie (in projection at least) 
within the SNR shell and are potential counterparts for 3EG 
J1410-6147. 

Given the presence of young, energetic pulsars, the SNR 
and a putative PWN, we undertook radio and X-ray obser¬ 
vations of this region. The radio observations were primarily 
aimed at determining an H l kinematic distance to the SNR. 
A known distance then enables us to consider the energet¬ 
ics of the system. In X-rays we hoped to detect the young 
pulsars and hence unravel the nature of the SNR and its 
PWN. 


2 OBSERVATIONS 
2.1 Radio Observations 

Radio observations of SNR G312.4—0.4 were made with the 
Australia Telescope Gompact Array (ATCA; Prater, Brooks 
& Whiteoak 1992). The ATCA is an east-west synthesis tele¬ 
scope located near Narrabri, NSW and consists of six 22 m 
antennas on a 6 km track. Five of the antennas are moveable 
along the track, enabling observations to be carried out in 
conhgurations of different baselines between 31 m and 6 km. 
The ATCA can observe at two different frequencies simulta¬ 
neously and is capable of recording all 4 Stokes parameters 
at each frequency. 

Our observations consist of two separate pointings and 
were obtained over 6 observing sessions, with four different 
array configurations, 6D, 1.5A, 750C and 375. Details are 
given in Table 2. Each observation was for 12 hours. Antenna 


Figure 1. MOST image of SNR G312.4-0.4 at 843 MHz 
(Whiteoak & Green 1996), overlaid with EGRET confidence con¬ 
tours at 68%, 95% and 99% (Hartman et al. 1999). 

1 was not functioning during the observations on 2001 Jan 
9; hence, only 10 baselines were recorded. 

Our observations comprise continuum data split into 32 
channels, each of 4 MHz, for a total bandwidth of 128 MHz 
centred at 1384 MHz, and Hi data with a 4 MHz band¬ 
width and 1024 channels, centred at 1420 MHz. All Stokes 
parameters were recorded for the continuum data, enabling 
polarisation parameters to be calculated. The Hi observa¬ 
tions were made only in total intensity. 

2.2 X-ray Observations 

We obtained a 9.7 ks exposure of G312.4-0.4 with the Chan¬ 
dra X-ray Observatory on 2001 July 26, with pointing centre 
RA (J2000) 14 12 6.52, Dec (J2000) -61 45 44.8. The four 
ACIS-I chips, together with two chips from the ACIS-S ar¬ 
ray provide significant coverage of the SNR (see Fig. 5). The 
data were reduced with standard procedures from the GIAO 
2.2.1 package. No large background flares were observed dur¬ 
ing this exposure, so the full time interval was incorporated 
into the final maps. 
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2.3 Other Observations 

We have used several other data sets in this study: 

• Radio data from the MOST at 843 MHz (Green et al. 
1999). 

• Single dish continuum and H l data for this region, at 
21 cm, which were taken as part of the Southern Galactic 
Plane Survey (McGlure-Grifhths et al. 2001). 

• Images from the Gohen & Green (2001) study of this 
region of the Galactic Plane in the mid-infrared using obser¬ 
vations made with the Midcourse Space Experiment (MSX). 


3 RADIO DATA REDUCTION 

Reduction and analysis of the radio data were carried out 
with the MIRIAD package using standard techniques de¬ 
scribed in The Miriad User’s Guide (Sault & Killeen 1998). 
The primary calibrator PKS B1934—638 was used for flux 
density and bandpass calibration. The secondary calibrator 
B1329—665 was used to solve for antenna gains, phases and 
polarisation leakage terms. After calibration was completed, 
the continuum data was recorded as 13 channels each with 8 
MHz bandwidth. The spectral line data have 512 channels, 
each with a velocity resolution of 1.6 km s“^. 


3.1 Continuum Image 

The continuum data was reduced using mosaicing techniques 
in order to account for the two different pointing centres. 
The image was formed using a weighting scheme intermedi¬ 
ate between natural weighting (which minimises noise) and 
uniform weighting (which minimises sidelobe levels). This 
was achieved by setting the robust parameter to 0.5 in the 
invert task, giving nearly the same sensitivity as natural 
weighting but with signihcantly lower sidelobes. The image 
was deconvolved using a maximum entropy algorithm which 
is implemented in the task mosmem (Sault, Staveley-Smith & 
Brouw 1996). The final image has 3” pixels and a beamsize of 
25”. The rms noise in the image is 0.5 mjy beam“^, which is 
higher than the theoretical value but not unexpected consid¬ 
ering the excess background emission in the Galactic Plane. 
This image is shown in Figure 2. 

A second image was also formed, combining the ATGA 
interferometric data with single dish data from the Parkes 
radio telescope (McGlure-Grifhths et al. 2001). There are 
two possible methods for doing this. One is a non-linear 
technique whereby the data are combined in the uv plane 
and then jointly deconvolved. An alternate method, which 
is used here, combines the data in the Fourier domain af¬ 
ter both have been separately deconvolved (e.g. McGlure- 
Grifhths et al. 2001). This was implemented by the task 
iiranerge in MIRIAD. Including the single dish data adds 
the smoothly varying component of the source and the dif¬ 
fuse Galactic emission. 

Finally, a high resolution image was produced, with 6" 
resolution, using only the baselines with the 6 km antenna. 
This resolves out most of the extended emission and allows 
better sensitivity for point sources hidden beneath the ex¬ 
tended structure. The rms noise in this image is 0.07 mJy 
beam“^. 


3.2 Polarisation Images 

Stokes Q, U and V images were formed by inverting and 
cleaning each of the 13 frequency channels separately, us¬ 
ing the task pmosmem. A linear polarisation image, L, was 
formed for ea ch chann el from the hnal Q and U images, 
where L — . The position angle image, PA, was 

also calculated, where PA = itan“^(t7/Q). The L images 
were then summed over all channels to produce a hnal lin¬ 
ear polarization image and a bias was removed to ensure the 
noise in the image had zero mean. The hnal rms noise of this 
image is 0.14 mJy beam“^, signihcantly better than for the 
total intensity. The resulting polarisation image was divided 
by the total intensity image and clipped where the total in¬ 
tensity emission was less than 5cr (2.5 mJy beam“^), to form 
a fractional polarisation image. This is shown in Figure 3. 

3.3 Spectral Index 

In order to compute spectral indices using data from only 
two wavelengths, it is essential that the hux densities be 
accurately determined. For an interferometer, the integrated 
hux density measured depends on the length of the shortest 
baseline used, in other words, the largest spatial scale that 
is detected. It is, therefore, important that the data contain 
the same spatial scales (e.g. Gaensler et al. 1999; Katz-stone 
et al. 2000) . 

To match the ATCA and MOST data, we created a new 
continuum image at 20 cm, using only the pointing centred 
at RA (J2000) 14:12:00, Dec (J2000) -61:47:30. These data 
contain baselines < 750 m and have a similar beam (^ 40”) 
to the MOST image at 843 MHz. The two images still have 
diherent uv coverage, so spatial hltering was carried out to 
equalize the spatial content of the two images. Extended 
structure was excluded from the 843 MHz image and the 
1384 MHz image was smoothed to the lower resolution of 
the 843 MHz image. 

We used two methods for calculating spectral indices, 
spectral tomography and temperature-temperature (T-T) 
plots and the results were compared, for consistency. Spec¬ 
tral tomography produces maps of the following quantity: 

/t(at) = /i - (1) 

where at is the spectral index, as dehned by Si, oc i'°‘, and 
7i and I 2 are the intensities of a given pixel at v\_ (843 MHz) 
and V 2 (1384 MHz) respectively (Katz-Stone et al. 2000). A 
series of maps was made for —1.0 < at < 0.1 in steps of 0.1. 
Any regions of the image which vanish into the background 
have It{at) ~ 0 and a spectral index of at- 

T-T plots traditionally involve plotting the brightness 
temperature of two images against each other, pixel-by¬ 
pixel. However, the flux density in units of Jy beam“^ can 
be used providing the two images have the same beam size 
(Bock 1997). This method is independent of any zero error, 
and the slope of the line of best fit is given by 



where a is the spectral index. The uncertainty in the linear 
fit may be underestimated as adjacent pixels are not inde¬ 
pendent and there are typically many pixels per beam. This 
effect can be corrected for by rebinning each image, taking 
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Date 

Array 

config. 

Pointing centre (J2000) 
RA Dec 

(MHz) 

D2 

(MHz) 

1998 Nov 1 

6D 

14:13:22 

-61:33:31 

1384 


1999 Nov 27 

375 

14:13:22 

-61:33:31 

1384 


1999 Dec 17 

1.5A 

14:13:22 

-61:33:31 

1384 


2001 Jan 9 

750G 

14:12:00 

-61:47:30 

1384 

1420 

2001 Jan 12 

750G 

14:12:00 

-61:47:30 

1384 

1420 

2001 Feb 25 

375 

14:12:00 

-61:47:30 

1384 

1420 


Table 2. ATCA Observations of SNR G312.4—0.4. 


Figure 2. ATCA greyscale images of SNR G312.4—0.4 at 1384 MHz. The images have a resolution of 25”. The greyscale range of image 
(a) is 0 to 20 mjy beam“^ and that of image (b) is 0 to 30 mjy beam“^. The FWHM of the restoring beam is shown in the lower corner 
of each image. The positions of PSRs J1413—6141 and J1412—6145 are marked with white crosses and seven unresolved sources have 
been labelled as objects of interest as possible counterparts to the EGRET source. The western region of the SNR, part of which may 
be a PWN, has been enlarged to show more detail. 


every ith pixel where there are i pixels per beam in each 
dimension (Gaensler et al. 1999). T-T plots were made for 
subsets of the images at the northern rim of the SNR and 
also at the brightest part of the western region. A linear 
least-squares-fit was then used to determine the slope for 
each region. The results were checked for self consistency by 
varying the box size for the plots (i.e. the number of pixels) 
and by shifting the box by a few pixels. The slope stayed 
the same for each respective region and the fit improved 
with the inclusion of more pixels. 

3.4 Spectral-line 

After flagging and calibration of the two Hi datasets, we 
performed continuum subtraction (task uvlin) by fitting a 
second order polynomial to the line-free channels. A ‘dirty’ 
cube was produced using the task invert, selecting only 
those baselines which were longer than 125 m and using 


uniform weighting. The output cube consists of 150 channels 
each 2 km s“^ wide starting at —120 km s“^. The cube was 
not cleaned but was convolved with a gaussian restoring 
beam of size 45 ". 

We made a continuum image in an identical fashion 
using the two datasets and the same parameters in invert. 
In this case we cleaned the image appropriately and restored 
it with the same gaussian beam of FWHM 45 . 

The rms noise as a function of channel number was com¬ 
puted using the routine imstat. The sensitivity varies over 
the bandpass depending on the temperature of the H l emis¬ 
sion in a given channel. Although the presence of variable 
sky structure in the cube will adversely affect the computed 
rms noise, this was the best way to assess the significance of 
any detected absorption feature. The data were then scaled 
by the total flux in the continuum image. This allows the 
fractional absorption depth to be computed and the signifi¬ 
cance of the feature estimated. 
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An emission profile was obtained using the Parkes single 
dish data, for a location just outside the SNR (McClure- 
GrifBths et al. 2001). 


4 RESULTS 

4.1 Continuum Image 

Radio continuum images of SNR G312.4—0.4 with a resolu¬ 
tion of 25” and an rms noise of 0.5 mjy beam“^ are shown 
in Figure 2. The SNR is seen as an incomplete shell, with 
a horseshoe-like appearance. Its brightest regions are the 
extended western component and the north rim. The arc¬ 
shaped blow-out feature which is visible in the MOST image 
to the south of the SNR can only be seen very faintly here, 
because of partial primary beam correction. The brightest 
part of the western region of the SNR has a flux density 
of ^30 mJy beam“^. It extends over a region of 14' by 7' 
and has an integrated flux density of ^5 Jy, with the con¬ 
tribution from the shell subtracted. The morphology of the 
western bright region consists of a complicated filamentary 
structure, which is seen for the first time. The brightest parts 
of the northern rim of the SNR have flux densities ~20 mJy 
beam“^. 

There are several bright, unresolved sources in the im¬ 
age. Seven of the sources (labelled in Figure 2) are now dis¬ 
cussed. Source 1 (G312.11-0.2) is non-thermal with a flux 
density of 600 mJy at 408 MHz (Caswell & Barnes 1985) and 
166 ±3 mJy at 1384 MHz. It is classified as extragalactic by 
Cohen & Green (2001). Source 2a (G312.36-0.04) has a flux 
density of 173 ± 6 mJy at 1384 MHz and is classified as a 
compact H ll region (Cohen & Green 2001) but the nature of 
source 2b (flux density 75 ± 2 mJy at 1384 MHz) is unclear 
although it is probably non-thermal. The flux densities at 
408 MHz and 5 GHz given by Caswell & Barnes (1985) are 
affected by blending between the two sources. Sources 3 to 
6 have 1384 MHz flux densities of 74 ±6, 59 ± 6, 35 ± 7 
and 17 ± 8 mJy respectively. These sources lie within the 
SNR shell and confusion with the underlying SNR emission 
contributes to the large error bars. Source 5 is seen in the in¬ 
frared and although not classified by Cohen & Green (2001) 
is thought to be a compact H ll region. Source 6 has an X-ray 
and an infra-red counterpart and is thus likely to be stellar. 

The two pulsars listed in Table 1 are not visible in Fig¬ 
ure 2 but their positions are shown as crosses. They are de¬ 
tected in the high resolution image and our positions agree, 
within the uncertainties, with those from the Parkes Multi¬ 
beam Survey (Manchester et al. 2001). We measured flux 
densities at 1384 MHz of 0.79 ± 0.07 mJy and 0.73 ± 0.07 
mJy for PSRs J1413—6141 and J1412—6145 respectively. 

The long filaments to the north-west of the SNR have a 
flux density of ^20 mJy beam“^ and are also visible in the 
infrared (see section 4.5). 

There are several weak background sources in the image 
which are almost certainly extragalactic. These are unlikely 
to be responsible for the y-ray emission as their flux density 
is a factor ^ 100 weaker than the known EGRET blazars 
(Mattox et al. 1997). 


4.2 Polarisation 

The fractional polarisation image is shown in Figure 3. The 
image has been masked where the total intensity is less than 
2.5 mJy beam“^. The strongest polarisation comes from a 
small area at the southern tip of the western region of the 
SNR, which is polarised at around 5—10%. It is offset by 
about 1' from a bright spot in the radio continuum. This is 
an unusual feature, whose origin is unclear. 

The absence of any significant polarisation at 1384 MHz 
is somewhat surprising, particularly as the western edge of 
the SNR is polarised at the 20% level at 8.5 GHz (Whiteoak 
1993). The lack of polarisation at 1384 MHz can possibly be 
explained by beam depolarisation, resulting from changes in 
the degree of foreground Faraday rotation on scales smaller 
than the synthesised beam. Regions of significant beam de¬ 
polarisation at this frequency have indeed been determined 
for nearby regions of the Galactic plane (Gaensler et al. 
2001 ). 


4.3 Spectral Indices 

T-T plots were made for the western part of the SNR shell 
and the northern rim and the results are shown in Figure 4. 
As described in Section 3.3, only independent pixels have 
been used, in order to avoid underestimating the uncertain¬ 
ties. The two regions have significantly different spectral in¬ 
dices. The slope of the linear fit for each region is related to 
the spectral index as expressed in equation 2. The slope of 
the distribution of points along the rim gives a = —0.7T0.1, 
which is typical of synchrotron radiation seen from shell 
SNRs. The distribution of points from the western region 
gives a = —0.19T0.01. This value is much flatter than the 
typical spectral indices for shell SNRs (—0.4 a ^ —0.7) 
but lies well within the typical range expected for a PWN 
(0.0 < a > -0.3). 

Spectral tomography indicates a spectral index of be¬ 
tween —0.6 and —0.8 for the northern rim and a value that 
is much flatter for the western region, consistent with the 
T-T result. 


4.4 X-rays 

There are no bright sources in the X-ray image. An appli¬ 
cation of wavdetect to the 0.2—8 keV data using a thresh¬ 
old SN = 2 detected only 18 sources. Several of these were 
identified with bright stars; none correlate with the bright 
radio point sources. Presumably, most of these sources are 
coronal emission from field stars or (for a few hard absorbed 
sources) emission from background AGN. Any emission from 
the G312.4—0.4 complex should also be highly absorbed. We 
therefore formed a 2—8 keV image and replaced the counts 
from the apertures of all point sources detected in the full- 
band image with Poisson counts from the local 2—8 keV 
background. The background counts/pixel were generally 
quite constant across the frontside illuminated chips. We 
subtracted the average counts/pixel from chips 10—13 and 
S2 and a separate background from the back-illuminated S3 
chip. No exposure correction was applied as the fluctuations 
at the chip boundaries obscured any potential diffuse emis¬ 
sion. After binning and adaptive smoothing, the resultant 
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Figure 3. Fractional polarisation for SNR G312.4—0.4 at 1384 
MHz. The greyscale ranges from 0% (white) to 7% (black). Con¬ 
tours show the radio continuum and the levels are 5, 10, 15, 20, 
25, 30, 40 mjy beam“^. 

0.08 I-^^^^^^^-1 



0 I-^^^^^^^-1 

0 0.02 0.04 0.06 0.08 

Flux at 1384 MHz (Jy/beam) 

Figure 4. T-T plots for two sub-regions of SNR G312.4—0.4. The 
diamonds show points from the northern rim whereas the crosses 
show points in the western region. The slope of the best fit line to 
the distribution from the rim is 1.4±0.14, giving a = —0.7±0.1 
and the slope of the best fit line to the points from the western 
region is 1.1±0.07, giving a = —0.19±0.01. 

ACTS image is plotted in Figure 5 with the 20 cm radio con¬ 
tours overlaid. A few regions of excess emission are seen, but 
at the west limb of the SNR only a slight excess in the diffuse 
hard count rate is visible. No signihcant diffuse emission is 
seen below 2 keV, indicative of high absorption. 

This diffuse X-ray emission from the western region is 
marginally detected. In the (unsmoothed) 2—10 keV point 


Figure 5. Smoothed ACIS-I-I-S2/S3 (2—8 keV) image of unre¬ 
solved emission in the direction of G312.4—0.4. The radio con¬ 
tours and pulsar positions (crosses) are shown. Point source sub¬ 
traction was carried out and the nature of the remaining source 
is unknown. Only the region with clear structure is significant. 
The edges of the image are outside the chip coverage. 


source-removed image, we find 105±24 excess counts in a 
6' X 6' region covering the brightest portion of the radio 
image. For typical {uHi/rie) values 10-30, the dispersion 
measure inferred Hi column density to PSR J1412-6145 is 
comparable to or greater than the total column of Galac¬ 
tic gas in this general direction, Nh = 2.2 x 10^^cm“^ 
(Dickey & Lockman 1990). Adopting this latter column den¬ 
sity and assuming a power law distribution with an index 
F = 1.5, which is generic for the hard component of young 
pulsars and their nebulae, we calculate a flux in X-rays of 
fx(2 - 8keV) = 3.2 ± 0.7 x 10“^®erg cm’^ s"^ This flux 
value is fairly insensitive to Nh unless it substantially ex¬ 
ceeds 10^^cm“^. The ASCA archival images also show some 
hard X-rays associated with this diffuse source. There is en¬ 
hanced emission in the range 4—10 keV which has a flux 
of ~ 4.5 X 10“^®erg cm“^ s“^. This is consistent with the 
Chandra results. 

We have also looked for unresolved sources correspond¬ 
ing to PSR J1412—6145 and PSR J1413—6141. No counts 
are detected at either position. PSR J1412—6145 is close 
to the optical axis of the ACIS-I and the excellent point 
spread function allows us to put an upper limit on the 
counts (2—8 keV), assuming an unabsorbed energy distri¬ 
bution with a power law index F = 1.5, which is equivalent 
to a flux of 6 X 10“^®erg cm“^ s“^ (Su). At the position of 
PSR J1413—6141, the degraded point spread function allows 
an upper limit of 1.8 x 10“^^erg cm“^ s“'^(3(j), insensitive 
to a moderate range of Nh values. 
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Figure 6. 8 fim mid-infrared image with G" pixels, overlaid with 
radio contours at 20cm. The greyscale ranges from 0.78 MJy sr”*^ 
to 10 MJy sr“^. Contour levels are 5, 10, 15, 20, 30, 50, 90 and 
120 mjy beam“^. 


4.5 Infrared 

Figure 6 shows radio contours overlaid on the mid-infrared 
8 /tm image from the MSX satellite (Cohen & Green 2001). 
A very striking feature of Figure 6 is the presence of thermal 
and radio continuum filaments running diagonally across the 
image. The bright filaments in the infrared are slightly offset 
from those in the radio. Cohen & Green (2001) have made 
an extensive comparison of this region in the radio and the 
infra-red. They discuss the nature of the filaments (and other 
sources) in detail. 

Although there is some diffuse infra-red emission 
throughout the SNR it is clear that there are no major H ll 
regions within its boundaries. In particular the western part 
of the shell, which has a flat spectral index in the radio, is 
not an obvious thermal source. As previously mentioned, the 
radio sources 5 and 6 have infrared counterparts and there 
is also a cluster of strong infrared sources close to the top 
rim of the SNR which are likely to be H ll regions or stars. 


4.6 HI Spectral Line Data 

We compare our resulting absorption spectra with the emis¬ 
sion profile in that direction, obtained from the Parkes data 
(McClure-Grifiiths et al. 2001). Figure 7 shows the emission 
profile and directly underneath it, the absorption spectra for 
Source 1 (as marked in the continuum image in Figure 2). 
The brightness temperature is scaled to match the value in 
the Kerr et al. (1986) Hi emission survey of this part of the 
Galaxy. The dashed line is the la envelope, plotted as a func¬ 
tion of velocity channel and clearly shows that the noise in 
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Figure 7. Emission (top) and absorption (middle) spectra in the 
direction of the point source PI. The dashed line in the middle 
panel shows the Icr envelope as a function of velocity channel. 
The bottom panel shows the conversion of velocity to distance at 
this galactic longitude. 


the line is about 4-5 times higher than the rms noise off-line. 
Below the spectra we show the rotation curve appropriate for 
this longitude. To compute the rotation curve we use an an¬ 
alytic expression given by the best fit model of the Galactic 
rotation curve from Fich, Blitz, & Stark (1989). We adopt 
the standard lAU parameters (Kerr & Lynden-Bell 1986) 
for the distance to the Galactic centre (Ro = 8.5 kpc) and 
the solar orbital velocity (Vb = 220 km s“^). The rotation 
curve of the Galaxy is the azimuthally smoothed average of 
its velocity field and thus implicitly assumes circular orbits 
around the Galactic centre. 

The emission spectrum is somewhat unusual in that it 
shows relatively weak emission at zero velocity. High bright¬ 
ness temperature gas is seen out to the tangent point at 
—55 km s“^. At positive velocities we see 100 K gas out to 
-1-40 km s“^ and two further weak emission features at -1-70 
and +95 km s”*^. One would not expect to see absorption 
against these two weak features. 

The absorption spectrum for Source 1 (Figure 7) dis¬ 
plays deep absorption at all negative velocities out to the 
tangent point. Absorption with a similar optical depth to 
the negative velocity features is also seen from +20 to 
+50 km s“^. Source 1 shows absorption against all the main 
emission features and it is almost certainly extragalactic, 
confirming the classification of Cohen & Green (2001). The 
absorption spectrum for Source 2a is very similar to that of 
Source 1 and it seems probable that it is also extragalactic. 
This is at odds with the compact H ll region classification of 
Cohen & Green (2001) but it is likely that their infrared to 
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Figure 8. Emission (top) and absorption (middle) spectra in 
the direction of the western part of the SNR. The dashed line in 
the middle panel shows the la envelope as a function of velocity 
channel. The bottom panel shows the conversion of velocity to 
distance at this galactic longitude. 


radio flux ratio (on which their classification is based) is too 
high because of confusion with the (foreground) filamentary 
structure. None of the other point sources were of sufficient 
flux density to show significant absorption in their spectra. 

Figure 8 shows the absorption spectrum in the direc¬ 
tion of the western region of the SNR. The spectrum ex¬ 
hibits deep absorption features at negative velocities out to 
the tangent point. Features with a smaller optical depth are 
also seen, at 20 km s“^ and 40 km s“^, with a significance 
of 4 (t. Given the deeper absorption features out to the tan¬ 
gent point, we can certainly impose a lower distance limit 
of 6 kpc. Any conclusions beyond this depend then on the 
interpretation of the features at 20 km s“^ and 40 km s“^. 
If these absorption features are real, a lower distance limit 
of 14 kpc is implied. 

The absorption spectrum towards the northern rim 
is very noisy and cannot provide any direct evidence of 
whether or not the northern rim and the western region 
are physically connected. 


5 DISCUSSION 

Given the uncertainties in the interpretation of the absorp¬ 
tion spectrum, in the following discussions we adopt a dis¬ 
tance of bde kpc to the SNR complex. 


5.1 Origin of the SNR 

As the main shell is quite symmetrical, we estimate its centre 
to be at RA (J2000) 14 13 00, Dec (J2000) -61 43 00. The 
radius of the shell is then 13', equivalent to a physical radius 
of 23d6 pc. Assuming an ambient density of n ~ 1 cm“®, the 
mass swept up by the shell is ~1200dg M©. This implies that 
the SNR is well into the Sedov-Taylor (adiabatic) phase of 
expansion as it has swept up more than 20 times the mass 
ejected (Dohm-Palmer & Jones 1996). 

We can relate the initial conditions of the SN to the 
radius of the expanding shock front by 

R{t) = 5.39 pc (3) 

where Usi is the kinetic energy of the explosion in units of 
10®^ erg, t is the time after the explosion in kyr and n is 
the density of the medium into which the SNR is expand¬ 
ing in particles cm“®. Assuming that the age of the SNR 
is the same as PSR J1412—6145 (50 kyr), we find the ratio 
i? 5 i/n ^ 0.6dg. Alternatively, the age of PSR J1413—6141 
(13 kyr), yields a ratio of Usi/n ^ 9dg. Given the uncertain¬ 
ties of the pulsar ages and the SNR distance, both these val¬ 
ues are well within expected bounds for the adopted distance 
d = 6 kpc. However, if d > 14 kpc then PSR J1413—6141 
may be ruled out as the progenitor star for the SN explosion, 
as Esi/n becomes unreasonably large. 

The pulsar velocities can be calculated by assuming, in 
turn, that each was the progenitor star for the SNR. This 
implies a transverse velocity of u = 218dg km s“^ for PSR 
J1412-6145 and v = 288de km s"^ for PSR J1413-6141. 
Again, both these values are entirely consistent with pulsar 
velocities in general (Cordes & Ghernoff 1998), even if d > 14 
kpc. 

5.2 A PWN around PSR J1413-6141? 

We might expect to see a radio PWN surrounding PSR 
J1413—6141, a young, highly energetic pulsar. Gaensler et 
al. (2000) discuss two types of PWN which can form, a static 
PWN which is confined by the gas pressure of the surround¬ 
ing interstellar medium and a bow-shock PWN confined by 
ram pressure from the motion of the pulsar through the 
medium. The condition for a static PWN is given by 

nU® < 4 X 10^E34/t^ (4) 

where V is the pulsar velocity in km s”*^, E 34 is the spin- 
down energy in units of 10®'^ergs“'^ and the other symbols 
are as defined earlier. Taking a lower limit on n of 0.003 
cm“®, the velocity of the pulsar would need to be less than 
210 km s“^ for this condition to be met. If this is the case, 
the size of the static PWN can be computed via 

Rstatic = 0.14 pc (5) 

yielding a radius of 4.8 pc. This PWN would be resolved in 
our image for any sensible pulsar distance but distinguishing 
the PWN from the remaining SNR structure is virtually 
impossible. 

If the pulsar velocity is in excess of 210 km s“^, it 
will have overtaken its static PWN resulting in a bow-shock 
PWN with a characteristic scale size 
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0.63 



( 6 ) 


Even for a lower limit on n of 0.003 cm“®, the bow shock 
PWN would be unresolved in our image. Assuming a typical 
PWN spectral index a = —0.3 and integrating from 10 MHz 
to 100 GHz, if the PWN is unresolved then the efficiency of 
the pulsar is given by 


e = 4.8 X 10"® 


d^SiA 

E34 


( 7 ) 


where S' 1.4 is the PWN flux density at 1384 MHz in mjy. In 
our case we measure a flux density of 0.97 mJy at the pulsar 
position as compared to a catalogue value of 0.51 mJy for 
the pulsed flux. Generously assuming a PWN flux of 0.5 mJy 
the efficiency is then < 1.5 x 10“®dg. This is more than an 
order of magnitude less than expected for pulsars of this age 
(Gaensler et al. 2000). 

We therefore conclude that PSR J1413-6141 has a ve¬ 
locity less than 200 km s“^ and has formed a static PWN 
whose surface brightness is confused with other emission or 
below our detection limit. This again points to the fact that 
PSR J1413-6141 is likely not to be associated with the SNR. 


5.3 Evidence for a PWN near PSR J1412—6145? 

There is circumstantial evidence that the western part of the 
SNR is a PWN. Whiteoak (1993) reports 20% polarization 
at 8 GHz and, as we have derived, the spectral index of this 
region is rather flat. However, the lack of polarization at 
both 5 GHz and 1.4 GHz is puzzling. This either implies a 
substantial rotation measure or that the 8 GHz polarization 
is over-estimated. On the other hand, if this is a PWN we 
need to determine the nature of its power source. 

The integrated flux density of the PWN is ~ 5 Jy. Us¬ 
ing the measured spectral index of —0.2, we estimate the 
integrated radio luminosity from 100 MHz to 100 GHz to be 
^l.lxlO^'^dg erg s The PWN covers an area of 12 x 24 
dg pc and so we assume a value of ^9 pc for its radius in the 
subsequent discussion. We can relate the radio luminosity, 
Lu, of a PWN to the spin down energy of the central source 
via Lu = eE where e is the efficiency. Gaensler et al. (2000) 
showed that in young pulsars e was of order 10”“^, and hence 
one requires the central pulsar to have E ^ 10®®dg erg s“® 
in our case. A pulsar of this E should be highly visible in X- 
rays. Furthermore, the two known pulsars in the field have 
E ^ 10®® erg s“®, which is clearly not high enough to power 
a PWN of such size and luminosity. The lack of any obvious 
high E pulsar in the Ghandra image, and the low values of E 
for the known pulsars, are inconsistent with the calculated 
parameters for a PWN. If d ~ 14 kpc, then both Lu and 
the size of the region become improbably large for a PWN. 
The source of the electrons powering this radio enhancement 
remains unclear. 


5.4 X-ray emission 

The limits on the X-ray luminosity of PSRs J1412-6145 and 
PSRs J1413-6141 are 2.6x10®^ ergs“® and 7.8x10®® ergs“® 
respectively assuming a distance of 6 kpc. This implies an 
X-ray efficiency of less than 2.1x10“® and 1.4x10“®. These 
limits are similar to the expected values for these sources 


from the Possenti et al. (2002) relation and are also in line 
with other pulsars of similar ages, eg. the Vela pulsar and 
the recently discovered pulsar PSR J2229-I-6114 (Halpern 
et al. 2001). 

The X-ray luminosity of the diffuse emission near PSR 
J1412-6145 is 1.4x10®® erg s“®. This is comparable to the 
luminosity of the diffuse nebula surrounding PSR B1951-I-32 
in GTB 80, although the E of this pulsar is higher than that 
of PSR J1412-6145 by a factor of ^30 (Safi-Harb, Ogelman 
& Finley 1995). Therefore the weak diffuse X-ray emission 
seen from this region in the ASGA and GXO data could 
plausibly be a PSR J1412-6145 powered PWN. 


5.5 7 -ray emission 

The 3rd EGRET catalogue (Hartman et al. 1999) gives 
the integrated flux of the 7 -ray source as F = 64.2 x 
10“® photons cm“^ s“® with a photon index of 2.12. Using 
an upper cut-off of 10 GeV, this gives a 7 -ray luminosity, 
Lj = 2.1 X 10®®dg erg s“®, assuming isotropic radiation or 
1.7 X 10®®dg erg s“® into 1 steradian. Tompkins (1999) has 
computed a variability index for the EGRET sources, find¬ 
ing a value of 0.331®'for 3EG J1410—6147. Nolan et al. 
(in prep.) using a revised statistic similarly hnd the value 
0.321 q' 24. These results show only that the source is signif¬ 
icantly less variable than the bright EGRET blazars, sug¬ 
gesting a pulsar or diffuse emission source of 7 -ray flux. 

The latter value of L^ — 1.7 x 10®®dg erg s“® is 139% 
of the E for PSR J1412—6145 and 30% of the E for PSR 
J1413—6141, even at d = 6 kpc. Generally one observes an 
increasing efficiency with decreasing E and PSR B1055—52 
is the most efficient 7 -ray pulsar with efficiency of 15% 
(Thompson 1998) and E a factor of 10 lower than PSR 
J1413—6141. For PSR J1412—6145 to be the 7 -ray source, 
a small beaming angle and a high efficiency would need to 
be invoked. The efficiency of PSR J1413—6141 is an order 
of magnitude higher than expected for its E and therefore it 
is unlikely that either pulsar powers the 7 -ray source unless 
the distance is substantially smaller that that inferred from 
dispersion measure and Hi measurements. 

With the existence of a PWN uncertain and prima fa¬ 
cie evidence for 7 -ray emission from the pulsars unconvinc¬ 
ing, the origin of the 7 -ray source 3EG J1410—6147 re¬ 
mains a mystery. It is tempting to speculate on the exis¬ 
tence of a third young pulsar in this region. It would have 
E ~ 10®®'ergs“®, sufficient to explain the 7 -rays and power 
the radio PWN. The main problem with this explanation 
is that its X-ray efficiency would have to be extremely low 
(< 10“®) to explain its non-detection in our X-ray observa¬ 
tions. It would then then have the lowest value, by an order 
of magnitude, of Lx/L-y of any of the known pulsars. A plau¬ 
sible alternative model could be based on the interaction of 
the western edge of the SNR with a nearby molecular cloud 
(eg. Aharonian, Drury & Volk (1994) and Esposito et al. 
(1996)). 


6 CONCLUSIONS 

We have made a continuum image of the SNR G312.4-0.4 
and surrounding regions at 1384 MHz, with a resolution of 
25”. We also made observations in the 21 cm Hi line and 
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produced absorption spectra to obtain a lower distance limit 
to the SNR of 6 kpc. If absorption features at positive veloc¬ 
ities are real, the SNR is located beyond 14 kpc. Although 
the bright, western portion of the SNR resembles a PWN, at 
these distances PSR J1412-6145 seems inadequate to power 
the observed enhancement to the radio flux. The lack of a 
radio PWN around the 13 kyr PSR J1413—6141 indicates 
it is likely to be slow moving and hence not formed with 
SNR G312.4-0.4. We suggest that the 50 kyr pulsar, PSR 
J1412—6145, originated in the same explosion as the SNR. 

The two young pulsars superimposed on the SNR are 
not detected in the X-rays, although the implied 10“^ lim¬ 
its on the X-ray efficiencies are not strongly constraining. 
Neither of these two pulsars can account for the y-ray lu¬ 
minosity unless abnormally high efficiency, small beaming 
angle and a low distance are invoked. 

The most likely interpretation of the multi-wavelength 
data of this region is (i) the SNR is large, old (~ 50,000 yrs) 
and distant, with PSR J1412—6145 as the likely progenitor, 
(ii) the bright western region of the SNR is most probably 
an amorphous, flat spectrum part of the SN shell and (iii) 
the 7 -rays are unrelated to the SNR. 

The 7 -ray source 3EG J1410—6147 still remains uniden¬ 
tified. The launch of the Gamma-ray Large Area Space Tele¬ 
cope in 2005 will reduce the size of the positional error box 
to 1' and provide a resolution to this puzzle. 
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